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Abstract 

We show that for recently discovered large values of ^13, a superbeam with an average neutrino 
energy of ~ 5 GeV, such as those being proposed at CERN, if pointing to Super-Kamiokande 
{L ~ 8770 km) , could reveal the neutrino mass hierarchy at 5 cr in less than two years irrespective 
of the true hierarchy and CP phase. The measurement relies on the near resonant matter effect in 
the Vfj, — )• Ue oscillation channel, and can be done counting the total number of appearance events 
with just a neutrino beam. 
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I. INTRODUCTION 



The recent determination of ^13 by T2K [T] , Daya Bay [5] and RENO [3] indicates a value 
very close to the previous Chooz bound [4J. The results of Daya Bay/RENO 

sin'2^i3|n..H.. = 0.092(16)(5), 



I DayBay 

sin2 2^i3|j,ENO = 0.113(13)(19), (1) 

are in perfect agreement. These new results have now been incorporated in global fits of 
all available neutrino oscillation data [5, 6J, which indicate towards a non-zero value of 613 
at around 8 a C.L. and hint that sin^0i3 lies between 0.019 and 0.033 at 2cr CL. with a 
best-fit value of 0.026 for normal mass ordering [5j. Such large value of the 1-3 mixing angle 
opens up the possibility to discover leptonic CP violation and determine the neutrino mass 
hierarchy (MH), normal (NH) if Am|]^ = mg — ml > 0, or inverted (IH) if Am|^ < 0, with 
significantly less effort than previously thought. 

Many studies have been carried out in the last decade to establish the best strategy to 
perform these measurements p'J . CP violation must be searched for by comparing neutrino 
and anti-neutrino appearance oscillation probabilities in the atmospheric range {E/L ~ 
lAmgj^l). This can be realistically achieved using neutrino beams in the GeV range produced 
in accelerators. On the other hand the MH determination is a discrete measurement that 
relies on the MSW resonance effect [8l |9] , which is relevant in several accessible neutrino 
beams: atmospheric neutrinos [TU], supernova neutrinos [TT], and also accelerator neutrino 
beams for sufficiently long baselines [12] . Cosmology can also weigh neutrinos with precision 
and future CMB and LSS cosmological measurements have a chance to determine the light 
neutrino spectrum[T3l [H]. Finally, the observation of neutrinoless double beta decay with 
the next generation of experiments would not only imply that neutrinos are Majorana, but 
also that the hierarchy is inverted [15j. 

The possibility to measure the hierarchy with atmospheric neutrinos in the proposed 
iron calorimeter detector at INO [16j, combined with the results from accelerator beam 
experiments T2K and NOvA, has been revised recently in light of the large value of ^13 [T7] . 
A 3cr C.L. determination of the hierarchy is possible irrespective of the value of the CP 
phase, only for the most optimistic assumptions about a 100 kt INO detector and a running 
time of ~ 10 years. 
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Out of all these possibilities, the method that is probably less affected by systematic 
errors, or correlations with other unknown parameters, is the one that uses a neutrino beam 
produced in an accelerator. In particular, given the large value of ^13, the use of more intense 
conventional beams, the so-called superbeams, would suffice. The optimization of energy 
and baseline for the measurement of CP violation and the determination of the hierarchy are 
somewhat in conflict. The compromise often implies some destructive interference of both 
measurements: not knowing the hierarchy influences the measurement of the CP phase and 
viceversa fTE\ . 

A very good precision on the CP violating phase can be achieved by measuring the neu- 
trino and anti-neutrino oscillation probabilities at energies around 1 GeV and at relatively 
short baselines of L ~ 300 km in a megaton-size water Cerenkov detector [19], such as 
proposed in the T2HK experiment [20]. For these baselines matter effects are small and 
therefore the hierarchy cannot be determined, which affects the precision achievable in the 
determination of the CP phase. The hierarchy on the other hand becomes almost a digital 
measurement at the point of maximal resonant conversion, which implies higher energies 
and baselines of -E ~ 6 GeV and L < lO'' km. In this case the measurement can be done 
with just a neutrino beam. 

Instead of trying to compromise both measurements, one could devise an optimal exper- 
iment for each. In this paper, we show that a setup that can easily determine the hierarchy 
could be achieved by focussing one of the proposed superbeams at CERN or Fermilab, to- 
wards the existing and well-tested Super-Kamiokande detector. This could arguably be an 
efficient way to determine the hierarchy, if NOvA misses it. However one should keep in 
mind that the baseline CERN-Kamioka is 8770 km and this requires a dip angle of ~ 42°, 
which could be a technical challenge. Such very long baselines have been considered before 
together with significant upgrades in the detector technologies, in the context of neutrino 
factories [2T1I22], beta-beams [231 [21] and superbeams [2S]. The most important result of 
this work is to point out that an existing detector as Super-Kamiokande could do the job. 
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FIG. 1: In the inset, neutrino energy where P^e is maximum as a function of L for NH. At the outset, the 
corresponding maximum oscillation probability versus L. The results are obtained using the approximation 
of eq. ([2]), i.e. neglecting the solar splitting. 

II. OPTIMAL E AND L FOR HIERARCHY DETERMINATION 

As is well-known, the sensitivity of neutrino oscillation probabilities to the neutrino MH 
relies mostly on the matter effects in neutrino propagation^. Neglecting the solar mass 
splitting, the — )■ v^. oscillation probability for a neutrino of energy E that crosses the 
Earth with a baseline L is given by 

P^e = sin^ sin^ 2^13 sin^ (^^|^) (2) 
using constant line-averaged Earth matter density, with 

Am^^ = y^(Am2^ cos 2^13 - Af + (Am|i sin 2^13)2, 

(3) 

and A = 2\/2GFn^{L)E. Gp is the Fermi constant and rig is the electron number density 
in the Earth which depends on L. Let us first consider the NH, that is Aml-^ > 0. In 

^ For a similar discussion see [5^1 [37] . 
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order to maximize the neutrino oscillation probability, it is not only necessary to choose the 
resonance energy which maximizes the effective angle in matter 

9 ;t , „ Amoicos26'i3 
sm'2e,s\E=E... = h i?res=— (4) 

but also we need to maximize the oscillatory term by choosing L = L^ax [26] : 

n,{L)L\, = — . (5) 

iLmax tan 2^13 ^ ^ 

Using the best-fit value of sin^20i3 = 0.101 from [5], we find Lmax ~ 10^ km and Ej-^s ~ 6.6 
GeV. It is very important to note that the maximization of the probability would not have 
been achieved in terrestrial distances had 613 turned out to be smaller. 

If we impose such resonance conditions, the oscillation probability for the IH, Amg^ < 0, 
would be very much suppressed since 

m:sm'2e,s\E-E ^ ■ (6) 

13\E-E,,, 1 + 3C0S2 2^13 ^ ^ 

Since also for the IH case, the oscillatory term would be at most 1, the oscillation probability 
would be widely different for NH and IH: 

P^e(NH)|s,,„L_„ - ^, 

PM.....^< ,^,^„J,g^^ ^0M4. (7) 

This analysis shows that if we have a neutrino beam with E = E^^s and L = L^i^^, the 
measurement is just a digital measurement. Of course in real life choosing Lmax and E^^s is 
difficult. 

The existing baselines such as CERN-Kamioka {L = 8770 km) and Fermilab-Kamioka 
(L = 9160 km) lie somewhat below Lmax- In the inset of Fig. [l| we show the neutrino 
energy for which the — )■ u^, oscillation probability (as given by eq. ([2])) is maximum for 
different choices of baseline. Here we consider NH. At the outset, we depict the corresponding 
maximum oscillation probability as a function of L. In going from Lmax to the CERN- 
Kamioka baseline, we see that the probability is reduced by a few per cent from its maximal 
value, and the optimal energy remains almost same. 

An additional advantage of decreasing the basehne is that we get closer to the so-called 
magic baseline |2T] . This is the baseline where solar oscillatory terms are strongly suppressed 
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Central (true) Values 


External 1 a error 


sin^ 6*13 (true) = 0.026 


a(sin2 ^13) = 13% 


Am|i(true) = 2.53 x lO'^ eV^ (NH) 


a(Amii) = 4% 


Am|i(true) = -2.40 x lO'^ eV^ (IH) 


a{Aml^) = 4% 


sin2 6*23 (true) = 0.49 (NH) 


(7(sin2^23) =8% 


sin^ 6*23 (true) = 0.53 (IH) 


a(sm2^23) = 8% 


Am|i(true) = 7.62 x 10~^ eV^ 


cj(Am|i) = 3% 


sin^ 6*12 (true) = 0.32 


a(sin2^i2) = 5% 


/9(true) = 1 


a{p) = 5% 



TABLE I: Best-fit values of oscillation parameters & their la estimated errors. In the last row, p is the 
Earth matter density, relative to the value given by the PREM profile. 



(and therefore the dependence on the CP phase): 

ML)L\^^^^.^ = L^agic ^ 7690 km, (8) 

and this choice of magic basehne is independent of the neutrino energy. Note that, we have 
neglected the solar sphtting in the previous analysis. This is quite a good approximation 
for the large value of 6*13 that we consider and because of the fact that we are close to this 
magic baseline. 

The existing basehnes CERN-Kamioka and Fermilab-Kamioka happen to be in the right 
ballpark, where the matter enhancement is still large while being close to the magic baseline, 
the sensitivity to the hierarchy will be all most independent of the CP phase. It makes 
therefore all the sense to ask the question how well the existing Super-Kamiokande can 
resolve the neutrino MH in conjunction with long baseline superbeam. 

For all the simulation results presented in this paper, we calculate the exact three gener- 
ation oscillation probability using the full realistic PREM [2H] profile for the Earth matter 
density. Unless stated otherwise, for all calculations, we use the central (true) values of the 
oscillation parameters given in the first column of Table [T] which are obtained from the global 
fit of world neutrino data |5]. We vary 6cp{tTue) in its allowed range of 0-2tt and study both 
the mass hierarchies. In all fits, we marginalize over all oscillation parameters, the Earth 
matter density by allowing all of these to vary freely in the fit. We impose simple Gaussian 
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Band: CP phase |0,2ll] Band: Maltcr Density [+5% to -5%] Band: CP phase [0,2jl] & Matter [+5% to -5%] 




Energy [GeV] Energy [GeV] Energy [GeV] 



FIG. 2: P^e as a function of E for the CERN-Kamioka baseline. We have taken sin^ 2^13 ~ 0.101 and 
for aU other oscillation parameters we assume the central values of Table [ij In left panel, the band portrays 
the effect of unknown 6cp- In middle panel, the band shows the effect of ± 5% uncertainty in the PREM 
profile on P^g. The combined effect of unknown 6cp and ± 5% uncertainty in the PREM profile is depicted 
as a band in right panel. 

priors on these parameters, with the corresponding 1 a errors as mentioned in the second 
column of Table |l] which are taken from [5J. Note that, in our study, we have imposed a 
prior on sin^ ^i3(true) with the la error of 13% based on the information from [5j, but its 
impact is marginal. The external information on the Earth matter density (p) is assumed 
to come from the study of the tomography of the Earth EO]- In the fit, we allow for 
a 5% uncertainty in the PREM profile and take it into account by inserting a prior and 
marginalizing over the density normalization. The CP phase 6cp is completely free in the 
marginalization. 

In Fig. [2| we show the full three-flavor oscillation probability — )■ z/g using the PREM [2H] 
density profile for the CERN-Kamioka baseline as a function of neutrino energy. We allow 
^cp to vary in its entire range of to 2it and the resultant probability is shown as a band 
in left panel of Fig. |2| with the thickness of the band reflecting the effect of 5cp on P^e- 
Since this baseline is close to the magic baseline, the effect of the CP phase is seen to be 
almost negligible. This figure is drawn assuming the benchmark values of the oscillation 
parameters given in Table [Tj We present the probability for both NH and IH. As expected 
the probability for the NH is a bit lower than 1/2 but still close to this maximal value. The 
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probability for NH is hugely enhanced for the neutrinos, while for the IH matter effects do 
not bring any significant change which is the key to distinguish between NH and IH. In 
middle panel of Fig. [2| the band shows the effect of ± 5% uncertainty in the PREM profile 
on P^e- The combined effect of unknown 6cp and ± 5% uncertainty in the PREM profile on 
Pfj^e is depicted as a band in right panel of Fig. |2] 

III. EXPERIMENTAL SETUP 

A. CERN based Superbeam 

Superbeams are conventional neutrino and anti-neutrino beams created from the decay of 
horn-focused pions. This superbeam technique is well established and understood. Future 
possibilities of having high-intensity proton sources at CERN and their relevance for high 
intensity neutrino beams were discussed in detail in Ref . [31] . In the context of the planned 
LHC upgrades, the CERN accelerator infrastructures hopefully will be able to supply high 
power beams. One of the possibilities discussed in [31.j is a new high power accelerator (HP- 
PS2), with a proton energy of 50 GeV and a beam power of 1.6 MW, resulting in 3 x 10^^ 
protons-on-target (pot) per year. A conceptual design is being developed in the context of 
the LAGUNA-LBNO design study [321 ES]- We will take this design as our reference setup, 
although obviously there are other alternatives. The expected neutrino fluxes from this 
machine have been estimated for the CERN-Pyhasalmi baseline by A. Longhin in [3ll 135] . 
We have simply scaled these fluxes to the longer L = 8770 km baseline, as L~^. There is 
however room for beam optimization, specially in what regards the neutrino energy, which 
could probably be increased. We will take as our reference an integrated luminosity of 
5 X 10^1 pot. 

B. The Super-Kamiokande detector 

Here in our study, we consider the existing and well understood Super-Kamiokande de- 
tector [361 ET] with a fiducial mass of 22.5 kt. In this work, we only rely on total signal and 
background event rates and no spectral information has been used. To estimate the total 
charged current (CC) signal event rates due to z/g appearance in the true neutrino energy 
window of 0.5 GeV to 10 GeV, we use the pre-cut efficiencies as given in Table 1 of [38] . 
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in the context of the T2KK proposal. These efficiencies are estimated based on the crite- 
ria that the events are fully contained inside the fiducial volume, have a single Cerenkov 
ring recognized as electron-like, and with no Michel electron present. It gives an average 
efficiency of 37% in the 5 GeV to 10 GeV true neutrino energy window where we have the 
maximum near resonant matter effect. For z/g appearance channel, we have considered three 
different types of background: intrinsic Ve contamination of the beam (Int) with the same 
efficiency as of signal, the number of muon events which will be misidentified as electron 
events (Mis-id), and neutral current (NC) events. In order to calculate the effective number 
of Mis-id and NC backgrounds in the true neutrino energy window of 0.5 GeV to 10 GeV, 
we take the corresponding efficiencies from Table 1 of [38] . Just to give an idea, in 5 GeV 
to 10 GeV true neutrino energy window, the NC background rejection efficiency is 90%. 

In [35, further cuts are applied to improve the signal-to-noise ratio for the appearance 
signal. These however have been optimized for the lower-energy JHF beam, so it would be 
necessary to adapt them for the neutrino beam considered in this work. For example, one 
could easily think of further improvements in the reduction of the dominant NC background 
from kinematical cuts {e.g. setting lower and upper threshold of reconstructed energy). The 
use of multi-ring events would also most probably provide valuable information. Since such 
optimizations require a detailed Monte Carlo study of the Super-Kamiokande detector, we 
take the conservative approach of using only the pre-cut efficiencies from [38] and leave the 
possibility of using optimized kinematical cuts for a more detailed study. 

We also include the information coming from the z/^ disappearance channel with the same 
signal efficiency that we have considered for electrons. For this channel, NC events are the 
only source of background because the 'wrong-sign' contamination is negligible. The NC 
background to the muon signal is expected to be significantly smaller because the single- 7r° 
contamination is negligible in this case. We have assumed 10 times more suppression in the 
NC background compared to the z/g appearance case. In any case, the disappearance signal 
gives a very small contribution towards the total MH sensitivity. 
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Channel 


CERN-Kamioka (8870 km) 


Signal 


Background 


CC-1 ring 


Int+Mis-id+NC = Total 


z/^ ^ V, (NH) 


40 


1+2+16=19 


l/e (IH) 


2 


1+3+16=20 


i/^ ^ z.^ (NH) 


84 


2 


i/^ ^ v^, (IH) 


89 


2 



TABLE II: First two rows show the total signal and background event rates in the v^. appearance channel 
taking sin^ 26'i3 ~ 0.101 and (5cp = 0°. Here 'Int' means intrinsic beam contamination, 'Mis-id' means 
niisidentified muon events and 'NC stands for neutral current. The last two rows depict the same for 
v^^ disappearance channel with only NC background. Results are shown for both the choices of the mass 
ordering with total exposure of 5 x lO'^^ pot. 



IV. RESULTS 



A. Event rates for CERN-Kamioka baseline 

In Table [Tl| we show the expected number of signal and background event rates for 
appearance and disappearance channels for both the choices of the hierarchies using the 



pre-cut efficiencies of the Super-Kamiokande detector as described in Sec. |III B[ To estimate 
these event rates, we consider sin^ 26'i3 = 0.101 and 5cp = 0°. For all other oscillation 
parameters, we have used the benchmark true values as given in Table [1} It is quite evident 
from Table |Tl] that there is a huge difference between the expected neutrino signal event 
rates for NH and IH in the appearance channel due to the presence of near resonant matter 
effect in the case of neutrino with NH. As far as the background is concerned, the main 
background contribution comes from NC events while the background coming from other 
two sources are quite small as can be seen from Table [Hj For disappearance channel, the 
difference between the signal event rates for NH and IH is not very large leading to a very 
small contribution towards the MH sensitivity. 

In Fig. [sj we show the expected number of events (signal + background) for the CERN- 
Kamioka baseline in the appearance channel as a function of sin^2^i3. Here we vary 5cp 
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Band: CP phase [0,2n] 




sin 2013 

FIG. 3: The expected number of appearance events (signal + background) as a function of sin^ 2^13 for 
the CERN-Kamioka baseUne. The result is shown for a total luminosity of 5 x 10^^ pot. The band portrays 
the impact of unknown S(jp . The vertical dashed brown line shows the present best-fit value of sin^ 2^13 of 
0.101. 

in its allowed range of 0-2-7? and the band reflects its impact. For sin^ 26'i3, we take the 
2cr allowed range as suggested in 0. We present the results for both the choices of the 
hierarchies taking a total luminosity of 5 x 10^^ pot. The difference in the number of events 
for NH and IH is quite large even for the value of sin^ 2^13 as small as 0.05. The lower part 
of the red band suggests that for NH, the number of events increases from 45 to 60 while 
we increase sin^ 2^13 from 0.05 to 0.14. 



B. Measurement of the Mass Ordering 



Here we explore the capability of the set-up described in Sec. HI to make a measurement 
of the true neutrino MH. A 'discovery' of the MH is defined as the ability to exclude any 
degenerate solution for the wrong (fit) hierarchy at 5 a confidence level. For our sensitivity 
calculation, we include a 5% systematic on the total number of signal events and a (uncor- 
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Normalized Exposure 



FIG. 4: Ax^ for the MH discovery as a function of the normahzed exposure. Here, 1 in x-axis corresponds 
to 5 X lO'^^ pot. The band portrays the effect of the unknown (5cp(true). We have assumed sin^ 26'i3(true) = 
0.101 and for all other oscillation parameters, the central values are taken from Table [l] 

related) 5% systematic on the total number of background events. They are included using 
the pull method as described in e.g. reference |39l HQ]. We perform the usual analysis 
using a Poissonian likelihood function adding the information coming from z/g appearance 
and disappearance channel. 

The total Ax^ obtained for the wrong hierarchy hypothesis is shown in Fig. |4] as a 
function of the exposure normalized to our reference exposure of 5 x 10^^ pot. The two 
bands correspond to the true NH and IH cases respectively, the width of the band in each 
case corresponds to the variation of the true value of the phase Sqp E [0, 27r]. Fig. |4] shows 
that to measure the neutrino mass ordering at 5 a C.L. irrespective of the choice of true 
hierarchy and the value of 6cp{tTue), we need an integrated luminosity of 5.4 x 10^^ pot. To 
achieve the same at 3(T C.L., a total luminosity of 1.73 x 10^^ pot is required. 

In Fig. [5} we address the issue that how the sensitivity of the proposed set-up towards MH 
discovery will change if we vary sin^ 2^13 (true) in its allowed 2 a range^. Here we consider 

^ We have used the same error on ^13 while changing its true value. This might not be the case in real 
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CERN-Kamioka Superbeam 
MH discovery (rate based analysis) 




sin 29 13 (true) 

FIG. 5: Ax^ for the MH discovery as a function of sin^ 26'i3(true). Here we consider a total exposure of 
5 X 10^^ pot. The band reflects the effect of the unknown (5cp(true). The vertical dashed black line shows 
the present best-fit value of sin^ 26*13 (true). 

our reference exposure of 5 x 10^^ pot. It is quite evident that with this exposure for true 
NH, this set-up can exclude IH well above 5 a CL. for any values of 5c7p(true) in the allowed 
2cr range of sin^ 2^13 (true). We can accomplish the same at 4.4 a CL. for true IH. 

Here we would like to point out that there is probably margin for optimizing the signal- 
to-noise ratio further. For example, we have checked that increasing the threshold energy 
from 0.5 GeV to 4GeV, improves the sensitivity. Considering NH as true hierarchy we need 
2.9 X 10^^ pot (3.4 X 10^^ pot) in the case of 4 GeV (0.5 GeV) threshold energy to exclude IH 
at 5(T C.L. irrespective of the choice of 5c'p(true). However, a detailed Monte Carlo study 
is needed to perform this optimization in reality. 

Other baselines could be realistically considered such as Fermilab-Kamioka (9160 km) 
assuming the same superbeam would be build at Fermilab. We find slightly worse but very 

scenario. 
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similar results for such longer baseline, everything else being the same. A more significant 
improvement is found for the shorter CERN-INO baseline (7360 km). In this case, however 
the detector would be an iron calorimeter, where it is more challenging to see electrons in the 
few GeV range. Reducing the thickness of iron plates in the ICAL@INO detector similar 
to the MINOS detector may create an opportunity to observe the electrons. A detailed 
detector simulation study would be needed to estimate the efficiencies and backgrounds in 
that case. 

V. CONCLUDING REMARKS 

In summary, we have explored a new strategy to determine the neutrino MH that employs 
a superbeam, as those being proposed at CERN, focused towards the existing and well- 
understood Super-Kamiokande detector, 8870 km away. We have shown that the neutrino 
beam resulting from a proton source with a total exposure of 5.4x10^^ pot can reveal the 
neutrino MH at 5 o" irrespective of the true hierarchy and CP phase, via a simple counting 
of z/e events at the far location. We believe there is margin for optimization both in the 
neutrino beam energy and in the detector signal-to-noise ratio. The proposed measurement 
of the neutrino MH would probably come earlier, and therefore fit nicely as a first step in 
the long-term plan to measure the more subtle effect of the CP violating phase. The latter 
requires a more ambitious program involving, for example, a superbeam running in both 
polarities, together with a significant upgrade in the detector technology, e.g. a factor of 20 
increase in the size of Super-Kamiokande, as proposed in the T2HK project [20j, or a more 
challenging detector technology such as C(20) kt Liquid Argon detector [IT]. 
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